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ABSTRACT 

Context. Although Blazars are a small fraction of the overall AGN population they are expected to be the dominant population of extragalactic 
sources in the hard X-ray and gamma-ray bands and have been shown to be the largest contaminant of CMB fluctuation maps. So far the number 
of known blazars is of the order of several hundreds, but the forthcoming AGILE, GLAST and Planck space observatories will detect several 
thousand of objects of this type. 

Aims. In preparation for these missions it is necessary to identify new samples of blazars to study their multi-frequency characteristics and 
statistical properties. 

Methods. We compiled a sample of objects with blazar-like properties via a cross-correlation between large radio (NVSS, ATCAPMN) and X- 
ray surveys (RASS) using the SDSS-DR4 and 2dF survey data to spectroscopically identify our candidates and test the validity of the selection 
method. 

Results. We present the Radio - Optical - X-ray catalog built at ASDC (ROXA), a list of 816 objects among which 510 are confirmed blazars. 
Only 19% of the candidates turned out to be certainly non-blazars demonstrating the high efficiency of our selection method. 
Conclusions. Our catalog includes 173 new blazar identifications, or about 10% of all presently known blazars. The relatively high flux 
threshold in the X-ray energy band (given by the RASS survey) preferentially selects objects with high /x/./i ratio leading to the discovery of 
new High Energy Peaked BL Lac (HBLs). Our catalog therefore includes many new potential targets for GeV-TeV observations. 
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1. Introduction 

Blazars are the rarest (~ 5%) and most extreme type of Active 
Galactive Nuclei (AGN) known. Historically, the classification 
of AGN has been largely based on observational characteristics 
leading to the proliferation of different classes of objects. All 
sources of this type, however, can be seen as part of a general 
paradigm in which AGN are divided into Thermal Emission 
Dominated (TED) AGN, where the emitted radiation is mostly 
generated through the accretion process onto a super-massive 
black hole, and Non-Thermal Emission Dominated (NTED) 
AGN, where the observed emission is mostly non-thermal and 
is generated in a jet of material moving away from the nucleus 
at relativistic speeds (e.g. |Giommi & Colafrancesco|[2006| l. 
Within this framework blazars are the small subset of NTED 
AGN in which the jet is closely aligned to the line of sight 
causing their emission to be strongly amplified by relativistic 
effects (as originally proposed by Blandford & Rees|1978 1. 

We recognize a source as a blazar if it shows the properties 
usually associated to aligned beamed emission such as strong 
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and rapidly variable emission in all energy bands, from radio 
to GeV, sometimes TeV energies, core dominated radio emis- 
sion with flat radio spectral index, superluminal motion of ra- 
dio compact regions, the presence of one sided jets (a jet on 
the other side is thought to exist, but with emission that is de- 
amplified by relativistic effects) and high brightness tempera- 
tures [Th ~ lO" - 10^^ K), close to or above the Compton limit 
(Tb ^ 1012). 

Sources that initially show only some of these properties, 
in later observations often also show the others, strengthening 
the hypothesis that these are basically equivalent and insep- 
arable features related to the same underlying physical process. 

Blazars include BL Lacertae objects (BL Lacs) where we 
observe a non-thermal optical continuum with no or very weak 
emission lines, and Flat Spectrum Radio Quasars (FSRQs) 
which exhibit both strong naiTow and broad emission lines. 
Furthermore, BL Lacs can be distinguished by the peak of the 
synchrotron emission in their Spectral Energy Distributions 
(SEDs). Objects with synchrotron peak at low energy (typi- 
cally in the Infra-Red) are generally found in radio surveys 
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and are called LBLs (Low Energy Peaked BL Lacs), whereas 
the rarer objects with synchrotron peak in the UV/X-ray band 
are called HBLs (High Energy Peaked BL Lacs) and are 
mainly selected in the X-ray band, where the maximum of 



their synchrotron power is emitted (e.g. Padovani & Giommi 
[1995). 

In this paper we discuss a large sample of candidate blazars 
that we have assembled using a multi-frequency selection tech- 
nique based on the NVSS ( jCondon et aL||1998j, ATC APMN 
(ATCA catalogue of compact PMN sources in Tasker (|2000[l 



RASS ( Voges et al.l999 ; Voges et al.|2000l l and GSC2 ( [Lasker 
[1995, McLean et al. 2000) catalogs. 

We have assessed the quality of the sample using a sub- 
sample of 816 objects for which data from Sloan Digital Sky 
Survey - Data Release 4 (SDSS-DR4; ,Adelman-McCarthy 



et al.|2()06i, 2dF Galaxy Redshift Survey (2dFGRS; jCoUess 



et al. 2001 1 and 2dF QSO Survey (2dFQSO; Shanks etal. 2000 1 



are available. 

We refer to our final list of 816 sources as the ROXA cata- 
log, which stands for Radio Optical X-ray ASDC catalog. 

We have built ROXA as a complement to the work of 
ISowards-Emmerd et al.] ( [2003j [2004 , 2005) as we share their 
goal of significantly enlarging the existing blazar catalogs 
in order to understand blazars statistical properties and to 
select interesting objects for upcoming high energy astronomy 
missions like AGILE, and GLAST. [Sowards-Emmerd et al.| 
( 2005 | l, with a single-band approach, selected flat-spectrum 
radio sources and then used optical telescopes for spectro- 
scopic follow-up observations, whereas in this work we use a 
multi-frequency approach that requires no additional telescope 
time. Moreover, the selection of radio sources results mostly in 
the discovery of LBL sources whereas our flux limits, which 
are mostly driven by the relatively shallow X-ray sensitivity 
of the RASS survey, favor objects such as HBL since they 
emit the maximum of their synchrotron emission near or 
within the X-ray band. This bias towards HBL sources is 
strengthened by the Ultra- Violet excess requirement in the 
SDSS spectroscopical target selection ( Richards et al.||2002| . 
The presence of HBLs allows ROXA to contribute with new 
targets for GeV/TeV observations. 



ROXA is also an augmentation of 'Collinge et al. (20051 
who searched for extragalatic quasi-featurless objects in the 
Second Data Release of Sloan Digital Sky Survey (SDSS-DR2; 



Abazajian et al. 2004) with low proper motion in USNO-B 



Catalog ( Monet et al.,.2003) to built a BL Lac candidate sam- 
ple. However, further information in other bands is required 
to confirm classifications therefore we used a multifrequency 
selection algorithm to build the candidate blazar sample and, 
although radio and X-ray selections could introduce biases, 
the high efliciency of the selection method allowed us to find 
many new blazars. 

The blazar selection algorithm is described in Section 2. 
In section 3 we focus on optical counterpart investigation and 
source classification criteria. Finally, we discuss our results 



and present our future projects in Section 4. 

Throughout the paper we assume a cosmology with param- 
eters in agreement with the results of Wilkinson Microwave 
Anisoti-opy Probe (WMAP; [Bennett et aL][2003|): Hq^IQ km 



s ' Mpc ', Qm=0.3, Qa^O.7 ( [Spergel et aLj2003t . 



2. The selection method 

For the definition of our sample we relied on the availability of 
large catalogs of astronomical objects combined with on-line 
services offering simple access to finding charts and magnitude 
estimates. 

The method is similar to that used for the DXRBS ( |Landt| 
et al. 2001 Padovani et al. 2006) and the Sedentary Survey 
( Giommi et al.|1999 2005 | l and consists in three steps: 



a first cross-correlation between radio and X-ray surveys 
(the NRAO VLA Sky Survey (NVSS; [Condon et al.|1998l ), 
ATCAPMN (ATCA catalogue of compact PMN sources in 
^ sker (200^1) and ROSA T All Sky Survey (RASS; [\^ges 



et al.|1999||Voges et a\.\20m) ). 



2. For each radio/X-ray match, optical magnitudes were re- 
trieved from the Guide Star Catalog (GSC2; I Lasker|1995" 
[McLean et al.|2000| ) 

3. For all radio/optical/X-ray matches we calculated the X-ray 
to optical (ffox ) and radio to optical (»„ ) spectral slopes 
and took only sources with and a^o values within the 



blazar area (for more detail see e.g. Perknan et al. 1998 
[Giommi et aT][T999j [Landt et al.[[2001[ [Padovani et al.[ 
[2006l l. 

2.1. The X-ray - radio cross correlation 

The NVSS catalog of radio sources includes nearly two million 



objects above a flux limit of 2.5 mJy at 1.4GHz (Condon et al. 
[TWSl ), the ATCAPMN catalog includes 7178 objects and it 
is reasonibly complete down to ~ 100 mJy ( [Tasker[ |2000| l 
whereas the RASS X-ray catalog is a list of 124.735 objects 



detected during ROSAT All Sky Survey ( [Voges et al.][T999 
Voges et al.|2000| ). 



As for the case of the Sedentary Survey (Giommi et al. 



[1999) and the DXRBS ( |Padovani et al.[|2006| l, we used the 



EXOSAT/BROWSE software to cross-correlate the NVSS and 
the RASS catalogs. The cross-correlation used a radius of 1 
arcmin and resulted in 16.596 matches. 

We then restricted our sample to sources with Galactic latitude 
I b\ > 20° finding 12.988 candidates. We kep t only those 
with A„ < 2.5cr„, where cr„ - '\Jcrl + crj is the total 
(radio-i-X-ray) positional uncertainty and A„ is the actual 
distance between the radio and X-ray positions. This resulted 
in 9663 entries. 

We estimated the percentage of accidental spatial coin- 
cidences of unrelated objects by shifting the coordinates of 
all the sources in one of the catalogs by a fixed amount and 
re-running the cross-correlation with the same parameters. 
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This procedure led to the statistical estimation that only 1 
% could be spurious associations. 



2.2. Blazar candidate selection 

For each RASS-NVSS match we obtained the magnitudes 
of the optical candidate counterparts from the GSC2 catalog. 
When no GSC2 source was present within the few arc-seconds 
positional uncertainty of the NVSS radio source, we assumed 
that the counterpart is fainter than Jmag=19.5, approximately 
the flux limit of GSC2 ( [McLean et al.|2000| l. In the very rai-e 
cases when more than one GSC2 object was present, the opti- 
cally brightest object was initially assumed to be the counter- 
part. 

We used EXOSAT/BROWSE software to get 5 GHz fluxes 
via a cross-correlation with several radio catalogs, such as 
NORTH6cm ( |Becker et al.]|T99T]l, GB6 ([G regory et al.'ll996|l 
and Parkes-MIT-NRAO (PMN; |Griffiflr& Wright 1993). For 
sources with no 5 GHz flux information we estimated it as 
Fschz - FiACHz ^ i^) """, where a™ is 0.25, the assumed 
average value for flat spectrum sources. 




Fig. 1. The cCqx - ffro distribution for the whole sample of 7.662 
blazar candidates 



We corrected X-ray and optical fluxes for Galactic absorp- 
tion and the redshift-dependent K-correction. For sources with 
no redshift information we assumed a typical redshift equal to 
the average value found in sources with the same X-ray to radio 
flux ratio (f„/fr ). 

To build our candidate blazar sample we have calculated 
the optical to X-ray (a^x ) and radio to optical (a,-a ) spectral 
slopes and selected objects with a^o > 0.2 where most of the 
blazars are located in the cKqx -ffro diagram (see e.g. Giommi 
|etal.|1999| for details). 

We relaxed the radio-loudness condition a^o > 0.2 to a^o > 
0.1 (where a^o and are the usual effective spectral indices 
defined between the rest-frame frequencies of 5 GHz, 5000 A 
and 1 keV) for objects which had no optical counterpart in the 
GSC2, since a lower limit on Jmag results in a lower limit on 

The selection process led to a sample of 7662 blazar can- 
didates. Fig.[T]shows the ai-o - ffox diagram and Fig.|2]plots all 
7662 sources in Galactic coordinates. The non-uniform space 
distribution in this last plot reflects the sensitivity of the RASS 
survey, which was higher around the ecliptic coordinates, and 
the higher radio flux limit (hence lower source density) of the 
ATCAPMN catalog at dechnations south of -40". 

3. Classification of optical counterparts using 
SDSS and 2dF data 

For all sources in our sample we used the SDSS on-line ser- 
vices (SkyServer-Radial search) to download SDSS measure- 
ments (right ascension, declination and u' g' r' i' z' magni- 
tudes) of each optical object lying within 6 arcseconds from 
each NVSS radio position. 

We found that 2.353 RASS-i-NVSS sources have SDSS- 
DR4 optical counterparts (including some multiple optical 




Fig. 2. The Galactic coordinates of the 7662 blazar candidates 
are plotted in AitofF projection 



counterparts, MOCs ). 

We then used the SkyServer DR4 Spectroscopic Query 
(proximity with radius 10 arcseconds) and the DAS online ser- 
vices to acquire spectra. We used previous downloaded SDSS 
positions to get spectral parameters and survey files for objects 
classified as galaxies, stars or "unknown" in the SDSS standard 
products. We found 669 spectra among the 2353 SDSS optical 
counterparts. 

In a similar way we searched for optical counterparts also 
in the 2dFGRS and in the 2dFQSO databases. This resulted in 
spectra and Bj magnitude estimates for additional 199 objects 
from 2dFGRS and for 94 from 2dFQSO (including MOCs). 
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3.1. Sources Classification 

All candidates that had MOCs were associated to a single op- 
tical object through a visual inspection procedure using the 
NVSS, ESO and NED online services. For each source we 
cross-checked the optical (ESO and SuperCOSMOS Surveys) 
finding chart with the RASS and NVSS positions and with the 
position in the FIRST catalog ( R.L. et al. 1997| l when available. 
In this way we were able to chose the best optical counterpart 
for each blazar candidate. Moreover, we excluded those objects 
with RASS source extent greater than 30 arcsec. We did not 
apply this condition to those detections with 20 or less photons 
since in these cases the uncertainty in the estimation of X-ray 
extent is too large. 

We built SEDs and analyzed optical spectra, with CaH&K 
evaluation when needed, to classify each object (the equivalent 
width of emission lines in radio-loud AGN depends strongly 
on the contribution from non-thermal jet emission, and, 
therefore, on the Ca H&K break value ( |Marcha et~aL]|1996| 
Scarpa & Falomo 1997 1). Sources with SDSS magnitude 
estimates, although not targeted for spectroscopy by SDSS 
pipelines, were kept in the sample only if already classified in 
the literature as reported in NED. 

SDSS and 2dF automated redshift estimation and spectral 
classification proved to be unreliable for several BL Lac 
objects because of the lack of emission lines. When possible, 
we re-calculated redshifts using the IRAF analysis package 
( |Tody|1986|[T993] l. 



We also defined a transition class for those objects 
that show properties in between two standard classifications. 
Regarding Radio Galaxies and BL Lacs we defined sources to 
be: 

- BL Lac if > 10"^"^ erg s ' or CaH&K < 0.4 or both 

- Radio Galaxy if CaH&K > 0.4 AND L, < 5 x 10"*^^ erg s ' 

- Radio Galaxy/BL Lac transition object if CaH&K > 0.4 

erg s" 



AND L, between 5 x lO''^ and lO'*'* 



Moreover we identified other transition classes for those 
objects showing optical properties which are borderline for 
their standard classification. These are BL Lac/FSRQ transi- 
tion object and R.G./FSRQ transition object. 

We re-calculated ofi-o and aox parameters estimating Vmag 
as in the following: 

- for SDSS objects: V = g' - 0.55{g' - r') - 0.03 

- for 2dF objects: V ^ Bj + 0.5 

- for object detected by both SDSS and 2dF we used SDSS 
magnitudes. 

Our final catalog ROXA includes 816 objects, 510 of which 
are confirmed blazars (i.e. 62.5%). Of the remaining sources 
1 10 (or 13% of the total) are confirmed QSOs (by the SDSS or 
2dF optical spectra) but their radio spectral slope is currently 
not available and therefore they remain blazar candidates. Only 
19% of the candidates turned out to be definitely non-blazars. 

Table [T] summarizes the statistics of our source classifica- 
tion. The complete catalog is reported in Table |2] The first 



column of Table [T| contains source classification, the second 
column gives the number of previously known objects (that 
is sources previously reported in the literature as blazars), the 
third column gives the number of new objects (classified in this 
work for the first time), the fourth column gives the total num- 
ber of objects and the last column gives the percentage with 
respect to the entire sample . 

The label confirmed blazars includes BL LAC, FSRQ, BL 
Lac/FSRQ transition object and BL Lac candidate. 

The label "others" is used to comprise several non-blazar 
sources, such as stars and variuos types of galaxies (Normal, 
Starburst and Seyfert). 

We point out that we identified 173 new blazars. 

Fig.[3]reports the multi-band indexes distribution of ROXA 
confirmed blazars. 




Fig. 3. ROXA a^o - Cox distribution of the confirmed blazars. 
Symbols are as follows: ★ - BL Lac, A - FSRQ, • - HFSRQ, e 
- BL Lac/FSRQ, -H - BL Lac candidate, x - HFSRQ candidate 



4. Conclusions and future work 

In this paper we presented the ROXA (Radio Optical X-ray 
ASDC) catalog, a list of of 816 objects including 510 con- 
firmed blazars. This catalog has been extracted from a sample 
of 7662 objects selected, through the cross matching procedure 
between the NVSS, ATCAPMN, RASS and the GSC2 catalog, 
to have properties similar to those of previously known blazars 
in the Oi-o -ffox diagram. 

We tested the selection method, also used in the DXRBS 
and in the Sedentary surveys, via a cross-correlation of the can- 
didate sample with the SDSS DR4 and 2dF surveys for which 
optical spectra are available. We analyzed SEDs and spectra to 
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Known objects 


new objects 


Total 


Sample % 


BL LAC 


182 


60 


242 


29.7% 


FSRQ 


147 


99 


246 


30.1% 


BL Lac/FSRQ transition obj. 


1 


7 


8 


1.0% 


BL LAC candidate 


7 


7 


14 


1.7% 


confirmed blazars 


337 


173 


510 


62.5% 


R.G./BL Lac transition obj. 





24 


24 


3.0% 


R.G./FSRQ transition obj. 





2 


2 


0.2% 


Radio Galaxies 


24 


8 


32 


3.9% 


SSRQ 


59 


43 


102 


12.5% 


QSORL 


60 


50 


110 


13.5% 


Galaxies NELG 


12 


2 


14 


1.7% 


BLRG 


2 





2 


0.2% 


others 


18 


2 


20 


2.5% 


Total 


512 


304 


816 


100% 



classify the candidate blazars and we confirmed the validity of 
the selection algorithm. 

Our final catalog consists of ~ 62.5% confirmed blazars, ~ 
13% QSO with no radio spectral information and of ~ 19% of 
sources that are surely non blazars. Since the ratio between flat 
and steep spectrum QSO is 254/102 (see tab. [T]) the majority 
of QSO with no radio spectral information is expected to be 
blazars. We can therefore conclude that our selection algorithm 
has an efficiency of (~ 70%). We also contributed several new 
blazar identifications finding 173 new blazars. 

Our method for the construction of the ROXA catalog is 
complementary to the work of Sowards-Emmerd et al. ( 2005[ l 
and Collinge et al. ( 2005| l as these authors used single-band 
approaches whereas we used a multi-frequency selection al- 



gorithm. Sowards-Emmerd et al. ( 2005 1 selected flat spectrum 
radio sources and then requested telescope time for follow- 



up spectroscopic observations whereas Collinge et al. (2005i 
searched for extragalatic quasi-featurless objects in SDSS-DR2 
with low proper motion in USNO-B Catalog to built a BL 
Lac candidate sample. However further information in other 
bands are required to confirm the classifications. Moreover, 
the selection of radio sources results in finding mostly LBLs 



by , Sowards-Emmerd et al. ( 2005j ) whereas our algorithm have 
greater chances of finding HBLs, because of our relatively high 
X-Ray flux limit that preferentially select objects in which the 
maximum of synchrotron power is emitted at high energy. 



We share the goal with Sowards-Emmerd et al. ( 2005 i to 
create as large as possible blazar catalogs to be used as targets 
for observations during the next planned new high energy 
missions (e.g. AGILE and GLAST). 

We plan to request high frequency radio observations in or- 
der to measure the spectral slope of the 110 QSOs with ra- 
dio measurements at a single frequency and to better evaluate 
the nuclear spectra of the 102 Steep Spectrum Radio Quasars 
(SSRQs) as their steep slopes could be due to extended radio 
components. 



In principle, by applying this selection method to the entire 
sky and not just to Umited regions, we could build a much larger 
sample suitable for the identification of foreground sources in 
CMB maps and in gamma-ray data. At present, the most suit- 
able optical data available are from the SDSS project. For this 
reason we intend to update the ROXA catalog as new SDSS 
releases become available. 
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